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We provide an extensive discussion on a scheme for Hamiltonian tomography of a spin- 
chain model that does not require state initialization [Phys. Rev. Lett. 102, 187203 
(2009)]. The method has spurred the attention of the physics community interested in 
indirect acquisition of information on the dynamics of quantum many-body systems and 
represents a genuine instance of a control-limited quantum protocol. 
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1. Introduction 

Recently, there has been a significant theoretical interest in the formulation of 
experimental-friendly criteria for the acquisition of informati on on the d etails of a 
given interaction model within a quantum many-body systemP^EIQISE] The moti- 
vations behind this are manifold. On one hand, the ability to manipulate and control 
quantum systems made out of a few elementary constituents (such as multi-photon 
states, multi-spin states in liquid nuclear magnetic resonance or cavity-implanted 
superconducting devices) has greatly improved in the last ten years, making the 
possibility to explore the dynamics of quantum few-body systems a reality. Clearly, 
the effects to test with the help of one of such systems would heavily depend on 
the arrangement of the correct and wanted form of interaction. On the other hand, 
many of the most advanced protocols for quantum information processing and com- 
munication rely on a fine pattern of coupling strengths across a lattice of interacting 
particles. This is important, for instance, when the task is to transfer a quantum 
state across a spin chainf^ It is thus very interesting to find a way to certify with 
a reasonable degree of confidence that the set of coupling strengths within a given 
experimental device is the one required by the associated theoretical protocol. 

From a practical viewpoint, it would be highly desirable to test if this is the case 
before running a protocol and to keep the invasiveness of such a test to minimum 
levels. A few important theoretical steps have been performed in this direction. In 
Ref. [U Burgarth et al. have proposed a scheme for the estimation of the paramo- 
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ters in a quantum spin-chain that is based on a clever application of the inverse 
eigenvalues problem, which is very well-known from classical inverse elastic theory.^ 
Along similar lines, in Ref. [T]we have proposed a scheme for "Hamiltonian tomog- 
raphy" , i. e. the determination of the coupling parameters in a chain of interacting 
spins, which requires only time-resolved measurements over a single particle, simple 
data post-processing and no initialization or prior knowledge on the state of the 
whole system. 

Our analysis takes advantage of an information flux approach^ that turns out to 
be particularly well-suited for the class of Hamiltonians to be probed by our tomo- 
graphic scheme. In fact, our method allows us to gather very important knowledge 
on the temporal behavior of the whole system simply by looking at the dynamics 
of a single element. From that, the parameters present in the Hamiltonian can be 
inferred. The protocol is efficient even when the spin-chain is affected by Markovian 
dissipativc and dephasing channels and can be generalized to Hamiltonian models 
that do not preserve the number of spin-excitations {i.e. interactions that do not 
commute with the total spin operator of the system). Very recently, Burgarth et 
al. have extended the basic idea of our Hamiltonian tomography method to general 
quadratic models, including the effects of transverse local magnetic fields and Ising- 
like tcrmsP^ Analogously, Wiesniak and Markiewicz have discussed Hamiltonian 
diagnostic tools that do not rely on state preparation.^ In this paper we provide 
a detailed account of the basic working principle behind the scheme discussed in 
Ref. [TJ The motivation behind Hamiltonian tomography, common to all the proto- 
cols put forward so far, is to provide a reliable way to fully characterize a many-body 
coupling model. As such, it can be seen as complementing and extending schemes 
that are already very well-known and successfully implemented, such as quantum 
process tomographyp] quantum state tomographjfl^ and detector tomographyff^ 
We hope that the non-demanding nature of the class of Hamiltonian tomography 
protocols suggested so far will soon spur the attention of the experimental commu- 
nity interested in many-body dynamics. 

The remained of the paper is organized as follows. In Sec. [2j we introduce the 
Hamiltonian tomography protocol for various classes of interactions. In Sec. [3j we 
use the scheme to determine the parameters in an engineered chain of eight spins. 
This serves as an illustration of the efficiency and working features of our protocol. 
Finally, we summarize our results in Sec. |4] 



2. Hamiltonian tomography for different spin-chain models 

The Hamiltonian tomography protocol proposed in Ref. [T] can be applied to vari- 
ous spin-chain models. A sketch of the scheme is presented in Fig. Q] The general 
scenario behind Hamiltonian tomography considers a chain of N spins mutually 
coupled according to a Hamiltonian model T-L. Of them, only the first spin, labelled 
hereafter as 1, is fully accessible and can be measured in any basis, the remaining 
sites being out of reach. Our task is to determine the coupling strengths entering T-L 
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Fig. 1. Sketch of the scheme for Hamiltonian tomography without state initialization, where M 
is the measurement performed on spin 1 and other spins are not accessible. 

by means of the information that we can extract from the measurements performed 
on 1 only. 

In Ref. [TJ we have discussed how the method can be used to infer the coupling 
parameters in JCY-like interactions described by the Hamiltonian 

AT-l 

Hi=J2 J i(XiX i+1 + YiYi+i), (1) 
i=i 

as well as the more general XY-like models 

JV-l 

n 2 = {Jx,iXiX i+1 + Jy&Yi+t). (2) 
1=1 

In the former (latter) case, J,: > {J a .i > 0, with a = X,Y) is the anti- 
ferromagnetic coupling strength of the pairwise interaction between qubits i and 
i + 1. The operators Xi, Yl and Z% denote the x, y and z-Pauli matrix of spin i, 
respectively. The condition of anti-ferromagnetic coupling can be dropped if the 
signs of all the parameters are known. It is also possible to extend the method to 
another class of Hamiltonians. In fact, using the results of Ref. 03J one can apply 
the scheme to determine the coupling strengths in the Ising-like model with a local 
transverse magnetic field given by 

JV-l N 
i=l i=\ 

Here, Jzs > is the anti-ferromagnetic coupling strength of the pairwise interaction 
between qubits i and i + 1 and Bi is the strength of the coupling of spin i to the 
local magnetic field affecting it. 

The key point of the scheme is that, for all the models described above, the time- 
dependent expectation values of simple single-spin operators contain information 
about the full set of coupling parameters. For instance, for the interaction model 
in Eq. ([T]), if the initial state of spin 1 is 

li^ = -^(10)^ ll^), (4) 
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we obtain 

(X 1 )(t)=±a 1 (t), (5) 

with a\(t) a function that depends on all the Ji's [see Eqs. (fTOj) and (fTTj) for a fully 
formal expression of such function]. On the other hand, for the model in Eq. ([2]), 
we need to consider the expectation values of two operators. The first one is again 
(Xi)(t), which provides the function a\(t) in a way completely analogous to what 
discussed above. As for the other one, we need to initialize spin 1 in an cigcnstatc 
of Yi 

li,}, = -^(10)^11)0, (6) 

so that we have 

(Y 1 )(t)=±0 1 (t). (7) 

Finally, for the Hamiltonian in Eq. we need to initialize spin 1 in one of the 
eigenstates of Z\, which we label {|0) 1 , |1) 1 }, so as to obtain 

(Z 1 )(t) = ±a 1 (t). (8) 

Notwithstanding the difference among the various models, the tomographic scheme 
required in order to infer the coupling parameters is the same. Such similarity 
holds in virtue of the fact that the evolution of the operators associated with spin 
1 depends on the operators of all the spins in the same way. This can be very 
effectively captured in a graphical way by using the information flux viewpoint put 
forward in Rcfs.[§]or by noticing that, adopting the language used in Ref. |T0J spin 
1 is infected by the rest of the chain in the same way, regardless of the details of 
the Hamiltonian model. 



3. Hamiltonian tomography at work 

In order to show the efficiency and working features of the scheme, here we ap- 
ply it to the case of a seven-coefficient problem. We have generated the set of 
randomly picked positive numbers {1.40,1.48,1.06,0.80,1.36,0.97,0.66}. Depend- 
ing on the model to be probed, they can be taken to correspond to specific 
parameters in the associated Hamiltonian. For instance, for the interaction in 
Eq. (Q]), they can simply correspond to the full set of J^'s. When Eq. ([2]) is con- 
sidered, they can be taken as the elements of a set of parameters constructed 
by taking Jx,i and Jy,i alternatively (for instance, we can identify them with 
{Jx.i, Jy,2, Jx.3, Jy,4, Jx,5, Jy,(n Jxj})- Finally, the random coefficients may cor- 
respond to the set of parameters {B\, Jz,i, -B2, Jz,2, -B3, Jz,3, B4} for a 4-spin chain 
evolving under the Hamiltonian in Eq. ([3]). 

For the sake of definiteness, we now describe the tomographic method in relation 
to the Hamiltonian in Eq. ([1]). It should be clear, though, that the corresponding 
results will hold for the other Hamiltonian settings, with due adjustments. Our 
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Fig. 2. Simulated dynamics of (Xi)(t) under the action of H with N = 8 and {Ji/J} = 
{1.40, 1.48, 1.06,0.80, 1.36,0.97,0.66}, sampled at steps of Jt = it /2b. The corresponding fit, per- 
formed using the trial function in Eq. is also shown. 

approach is as follows: we determine (X\)(t) at a few instants of time by simulating 
the evolution of the chain. This provides a discrete sampling of function a\ (t) which, 
as we said, contains full information on the coupling coefficients Jj's. In order to 
extract the values of such parameters, we shall fit the data points with a proper 
functional form. We have found that the trial function 

N/2 

a ( f ly) (t)=Y i A i cos(u H t) ( 9 ) 

i=l 

is in excellent agreement with the behavior of the simulated data. This allows us to 
determine the values of the amplitudes „4j's and frequencies u>i's. We now use the 
fact that ai (t) can be written as ^ 

*(*) = ( 10 ) 

1=0 

with S[ l) given by the recurrence formula 

5f = (-iy[j^ 1 sfz 1 1) + j j sf- 1 1) ], ai) 

where we have Jo = Jn = and the initial conditions 8^ = (1) for j ^ 1 (j = 1). 

Both Eqs. © and (TlTj|) are then expanded in Taylor series. By equating term by 
term the two series, one ends up with a linear system of algebraic equations whose 
solution gives the set of coupling parameters. Let us now illustrate the quantitative 
performance of the scheme: in Fig. [5J we show the dynamics of (Xi)(t) under the 
action of the Hamiltonian Hi constructed by means of the random set of coefficients 
given above. The expectation value (Xi)(t) is sampled at steps Jt = n/25 (the 
corresponding values are shown as dots in the plot). By fitting such data with 
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Eq. (|9]), we have obtained the following values for ,4j's and Wj's 

Ax = 0.0921, wi = 3.5929, 
■A 2 = 0.2748, cj 2 = 4.4941, 
A 3 = 0.3155, w 3 = 0.7821, 
Ai = 0.3176, w 4 = 1.6909. 

Having such TV values (that satisfy the constraint Yli=i = 1), we need TV — 1 
equations in order to find the TV — 1 parameters of the Hamiltonian under inves- 
tigation. The first TV — 1 non-zero terms in the Taylor expansion of a\ can be 
written as 

N/2 

^(-^(^E^- ( 13 ) 

We need to find the corresponding terms fj,j in the expansion of a±(t) as in Eq. (|10l) . 
Here we report only the first 4 terms (the other ones being rather cumbersome) 

Mi = - 2J X 2 , 

M2^(Jx 4 + ^Jf), 

4 (14) 
M3=-^J 1 2 [(J 1 2 + J 2 2 ) 2 + JlJ 3 2 ], 

M4 =^^{^i 6 + 3Jf J| + J^SJl + 2J|J|) + Jf[(jf + J|) 2 + J 3 2 J|]}. 
315 

By solving the system of equations r\j = /ij (j = 1,...,7), one finds the set of 
estimated coupling strengths 

j(eval) = 1 39998j j(eval) = 1 A8QQ ^ 
Aeval) = 1M0Q3 j(eval) = Q gQQQ5g 

J 5 (ei,oJ) = 1.36050, J 6 (et,ai) = 0.970524, 

j(eval) =0 660894j 

which are a very accurate estimation of the original ones. 



4. Remarks 

We have presented the working features of a scheme for Hamiltonian tomography 
that allows the identification of coupling parameters in various classes of spin-chain 
models through the study of the time dynamics of a single spin. The method is 
designed to work in a scenario of restricted access to the components of the chain. 
As no initial state preparation is necessary, one can perform measurements by in- 
terspersing them with the natural evolution of the system. Besides data acquisition, 
only a simple post-processing step is necessary: no conservation law associated with 
the interaction or a priori knowledge on the state of the system is required. Even 
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when a spin-chain is affected by environmental influences, the Hamiltonian tomog- 
raphy remains possible and reliable. Given the crucial role that proper coupling 
patterns have in the interference effects behind quantum many-body phenomena, 
non-demanding diagnostic methods, like the one proposed here, are important tools 
which need to be developed. 
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